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ABSTRACT: We develop and describe new approaches to the problem of interacting 
Fermions in spatial dimensions greater than one. These approaches are based on gener- 
alizations of powerful tools previously applied to problems in one spatial dimension. We 
begin with a review of one-dimensional interacting Fermions. We then introduce a simpli- 
fied model in two spatial dimensions to study the role that spin and perfect nesting play 
in destabilizing Fermion liquids. The complicated functional renormalization group equa- 
tions of the full problem are made tractable in our model by replacing the continuum of 
points that make up the closed Fermi line with four Fermi points. Despite this drastic ap- 
proximation, the model exhibits physically reasonable behavior both at half-filling (where 
instabilities occur) and away from half- filling (where a Luttinger liquid arises). Next we 
implement the Bosonization of higher dimensional Fermi surfaces introduced by Luther 
and advocated most recently by Haldane. Bosonization incorporates the phase space and 
small- angle scattering processes neglected in our model (but does not, as yet, address 
questions of stability). The charge sector is equivalent to an exactly solvable Gaussian 
quantum field theory; the spin sector, however, must be solved semiclassically. Using the 
Luther-Haldane approach we recover the collective mode equation of Fermi-liquid theory 
and in three dimensions reproduce the T^ In(^) contribution to the specific heat due to 
small angle scattering processes. We conclude with a discussion of our results and some 
speculation about future possibilities. 

PACS: 03.65.Fd, 11.40.Dw, 71.10.+X, 71.45.-d 



I. INTRODUCTION AND REVIEW OF A 1-D MODEL WITH 2-FOLD U(l) SYMMETRY 

Fermi liquid theory is now nearly forty years old. It is important to ascertain the range of 
its validity and determine whether more exotic generalizations, such as Luttinger liquids which 
exhibit spin-charge separation, exist. Shankar recently has emphasized the advantage of the 
renormalization-group (RG) approach over various types of mean-field approximations for answer- 
ing these questions^. In short, mean-field descriptions prejudice the outcome of the analysis by 
assuming that one, or at most a few, type of instability dominate the physics. RG analysis, on the 
other hand, treats all possible instabilities on an equal footing. Unfortunately, in spatial dimensions 
greater than one, the RG flows are described by nearly intractable functional equations. 

The approach we take in this paper is to slowly work up to the full problem by first reviewing 
rather well known one dimensional physics^. We then study a simplified model in two dimensions 
that incorporates some of the new physics that arises in higher spatial dimensions while still retain- 
ing the simplicity of one dimensional systems. Of course the price we pay for this simplicity is the 
drastic approximation to physical reality that we must make in order to arrive at the model: we 
completely neglect small angle scattering processes. Nevertheless, the model suggests a way to com- 
pletely reformulate Fermi liquid theory. Following Haldane's suggestion^ we now view Fermi and 
Luttinger liquids as zero temperature quantum critical fixed points characterized by infinite U(l) 
symmetry. The reformulation sheds light on how one might go beyond the drastic approximations 
of the model to include small angle scattering processes. 

Spin-charge separation occurs automatically in one spatial dimension, at least in the weak 
coupling limit and at long length scales. Consider the following low energy effective theory for 
excitations near the two Fermi points depicted in Figure [1]. The action in the non-interacting 
limit is given by: 

So = Jdx dt {^"a-V'La + ^]^d+ijRa} • (LI) 

Here L and R refer to the left and right Fermi points; d± = dt^ivfdx where Vf is the Fermi velocity 
which we will set equal to one for now. The Lorentz symmetry of this action guarantees that the 
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left moving Fermi fields are purely functions of the combination (x + ivft) whereas the right fields 
are functions of (x — iVft). The electron destruction fields Ca{x,t), where a =| or J, for up and 
down spins (with summation convention assumed), are related to the slowly varying continuum 
fields ipL,R by: 

c«(x, *) = 4f ie-*'^^" iJLaix, t) + e^'f^ i^R^ix, t)} . (1.2) 

Upon substituting this form into any given microscopic Hamiltonian the many-body interactions 
take the form of quartic and higher powers of the continuum fields. Most of the interactions are 
irrelevant in the renormalization group sense (see below) and the most general marginal interaction 
takes the form: 



+ Ac JlJr + As Jl^Jrr] 



1.3) 



Here for instance the charge current at the left point is defined by Jl ='■ tpitpLa '■ where the normal 
ordering symbols ":" indicate that we have subtracted the constant background charge density from 
the current to make < Jiixji) >= 0. The spin current is most conveniently expressed in matrix 
form: J£^(x) = ipj°'[x) ipLf3{x) — hS^ '>p'^^{x)ipL^{x). The matrix form can always be converted into 
the more familiar vector form with the identity: J£ = 2(^^)0 J lb where a = x, y, or z. Note that 
the spin current has no charge current component because it is traceless. It also has zero vacuum 
expectation value because in 1+1 dimensions the vacuum cannot break the continuous SU(2) spin 
rotational invariance by a quantum generalization of the Mermin- Wagner theorem. In the Fermion 
action, spin-charge separation is apparent even before Bosonization. That is, the interaction term 
involves only products of either pure spin or pure charge currents. The Gaussian part of the action, 
^o, also can be expressed purely in terms of separate products of the charge and spin currents (see 
below). 

Omitted from the action are terms that oscillate rapidly with wavevectors of order kf, interac- 
tions involving derivatives that arise from Taylor expansions of non-local interactions, and terms 
with more than four Fermion fields. Many of these terms break spin-charge separation; however, 
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each is irrelevant in the renormahzation group sense and the couphng constants flow rapidly to 
zero in the low-energy limit. To show the irrelevance, consider the scale transformation x ^ sx 
and t ^ st where s > 1. The Gaussian part of the action, ^o, remains invariant if we rescale 
the fields {ip^jip) -^ s~^''^{ip^ ,ijj). (Note that any non-linearities in the dispersion relation due to 
band structure are smoothed out as s ^ cxd.) Similarly, Sint remains invariant, showing that it is a 
marginal interaction. All other terms will, however, scale away at least as fast as an inverse power 
of s when s -^ cx). Thus, in one dimension, spin-charge separation occurs in the low-energy effective 
theory regardless of how the marginal interactions flow. Note that non-zero temperature acts as 
an infrared cutoff (since the time direction has a finite extent /? = j-^) that stops scaling towards 
the low-energy region beyond this scale. Irrelevant terms therefore persist at non-zero tempera- 
ture so the phenomenon of spin-charge separation must be construed as a zero-temperature critical 
property of the theory. 

Apart from the observation of spin-charge separation, the low-energy theory can be classified 
in terms of the symmetries that it obeys. In addition to the global SU(2) spin rotational symmetry, 
there exist two separate U(l) symmetries: one for each Fermi point. This U(1)l U(l)/j symmetry 
may be exhibited by considering the effect of separate left and right phase rotations by angles T^ 
and Fr on the Fermion variables: 

(1.4) 
tpRaix,t) -^ e*^« 1pR,a{x,t) . 

All of the currents are clearly invariant under this transformation, as the tp' fields transform with 
opposite phases. The physical meaning of the invariance is clear: the action, as it stands, conserves 
separately the number of left and right particles. We shall see that this special property has 
a natural generalization to higher spatial dimensions. Actually, one other marginal four-Fermi 
interaction can appear. The Umklapp term 

hiii^^Ri^Laf + H.c] (1.5) 

is permitted at half-filling in a periodic one-dimensional solid and it breaks the U{1)l tX" U{l)ji sym- 
metry down to the diagonal subgroup of ordinary U{1) transformations with F^ = F^j. It violates 
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the separate left and right U(l) symmetries because it transports two particles from one Fermi point 
to the other. Of course, total particle number remains conserved, and this conservation is reflected 
in the remaining diagonal U(l) symmetry. Like the other terms in the action, the Umklapp term 
preserves spin-charge separation because it transports charge, not spin, from one Fermi point to 
the other. To see this, note that it may be rewritten as: ^As {^af3WRVii){^'^^''pL^'4'L5) + H.c. where 
€a/3 is the totally antisymmetric tensor with ei2 = 1- Thus, only spin-singlet objects move from one 
Fermi point to the other. 

We now return to the problem without the Umklapp term and determine the RG flows and 
the nature of the fixed points. Bosonization of the Fermion fields is a powerful tool for addressing 
these questions. For now we use Abelian Bosonization^ and choose the spin quantization axis in 
the z direction. The current algebra will provide the vital link between the Fermion and Boson 
representations. We start by defining the normal-ordering operation carefully: 

JLaix,t) = : ipi{x,t) 7paix + e,t) : 

(1.6) 

= 1piix,t) lpaix + €,t)- < 1pi{x,t) ■lpaix + e,t) > . 

Here we place the spin index as a subscript on the V held to emphasize that we are no longer 
summing over it, and we imagine taking the e ^ limit at the end of our calculations. This 
"point-splitting" procedure regularizes ultraviolet divergences in our calculation. We choose real- 
space regularization because the connection between Bosons and Fermions occurs most naturally 
in real space. Momentum space regularization will be introduced later to permit the evaluation of 
momentum space integrals; differences between the two regularization procedures do not change the 
low-energy results. Currents for the right moving sector are obtained by making the replacement 
L ^ R and a simple calculation shows that left currents commute with right currents, while two 
left or two right currents at equal times obey the Kac-Moody algebra: 

[JLa{x) , JL/Siy)] = -^ <Jq/3 S' {x - y) . 

[jRa{x) , JRpiy)] = +Tr ^a/3 S' {x - y) . (1.7) 

ZTT 

6 



(To derive these relations, use the equal-time propagators for the Fermions < ^^"(a;)'(/'^^(0) >= 
dp 2^ ^-iid < ^Jf (x)'(/'_R/3(0) >= 6p 2^-) '^^^ coefficient of the derivative of the Dirac (5-function 
is known as the quantum anomaly. Note that it has the opposite sign for the left versus the 
right movers. The charge current defined previously may now be expressed in terms of these 
currents as: Jl(x) = JLlix) + JLi{x) and the z component of the spin current is simply: Jlz{x) = 
J^{x) — Jii{x). (However, the other two components of the spin current Jlx and Jiy are not so 
simply related to J^^ and Jl^ .) The charge and spin currents also obey the Kac-Moody algebra, 
but with twice the anomaly. 

We now introduce real- valued left and right moving free Boson fields (j)ia and (pj^ which satisfy 
the commutation relations: 

[(l>La{x) , (jiLpiv)] = -^ e(a: - y) , 

[(t)Ra{x) , (pRpiy)] = +1 e{x-y) , (1.8) 

and 

[(t^Ra{x) , (/>L/3(y)] = V4 , (1.9) 

where e(x) = 1 for a; > and = — 1 for x < 0. We also define canonical Boson currents 

and 

Jna{x,t) = +-j=^^^. (1.10) 

It is a remarkable fact that these Boson currents obey the same Kac-Moody algebra as the Fermion 
currents defined previously in Eq. [1.6]. To check the current commutation relation, take spatial 
derivatives of the Free Boson propagators: 

< </'La(x)0La(O) - (/'L(O) >= ^ In 



47r a + ix 



,2 -^ 1 , a 



< ^Ra{x)<t>Ra{0) " <pRaiO) >= y- In (1.11) 

47r a — IX 



(no sum over a in these expressions) to form the expectation value of the commutator. Here we 
see the first indication of equivalence between the Fermion and Boson representations. A further 
connection is revealed by examining the spectra of the free Fermion and Boson theories: both of 
which are linear and are in fact identical. To go further, we rewrite the quadratic Boson Hamiltonian 
in terms of the Boson currents defined in Eq. [1.10]: 



H = ^ Y. J dx ULix) + 4^ix)} . (1.12) 



The key point to be made here is that the currents appearing in this Hamiltonian could just as 
well be the Fermion currents Eq. [1.6] since the algebras are identical. Although in this represen- 
tation the Hamiltonian is a quartic function of the Fermion fields and as such would appear to be 
intractable, remarkably, as we have seen, it is equivalent to a free Fermion theory. 

We now consider the effect of the three types of interactions appearing in the action Eq. [1.3]. 
First, the current bilinears proportional to 5vc and 6vs simply renormalize the coefficients of the 
quadratic Boson Hamiltonian. [Although all three spin components of the spin current appear in 
Eq. [1.3], SU(2) invariance means that J1r{x)J^^{x) can be replaced by 3J|^(a;) so in this case it 
suffices to consider only the z-component of the spin current.] Using: 

we see that so far spin-charge separation is explicit as the renormalized Boson Hamiltonian may 
now be expressed as the sum of two pieces, H = H^ + Hg, that separately describe charge and spin 
excitations propagating at different velocities: 



and 



H, = -v, I dx Utix) + Jf,{x)} (1.14) 



Hs = ^vs Jdx UUx) + 4,(x)} (1.15) 

where Vc = 1 + Svc and Vg = 1 + Svg- It should be emphasized that the coefficients Svc and 5vs 
do not flow under RG transformations but instead simply renormalize the bare velocities. It may 



seem strange to have two different "velocities of liglit" - indeed, Lorentz invariance is broken. But 
since the charge and spin sectors are separate, Lorentz invariance is now manifest separately in the 
two sectors. 

Similarly, the charge current coupling, Ac, in the action Eq. [1.3] remains fixed as s ^ oo. To 
see this, note that it may be incorporated by adding another term quadratic in the Boson currents 
to the charge Hamiltonian: 

Hc= fdx {| [Jlix] + 4(x)] + Ac Jl{x) Jr{x)} . (1.16) 

To determine the effect of Ac on the spectrum, we must diagonalize this Hamiltonian. Diagonal- 
ization is accomplished via a Bogoliubov transformation that respects the Kac-Moody algebra. 
Let: 

Jl{x) = cosh(r/) J'i{x) + sinh(r/) J'^{x) 

and 

Jnix) = sinh(r/) J^(x) + cosh(r/) J'fi{x). (1-17) 

Then the primed charge currents obey the same algebra as the original (unprimed) charge currents. 
Upon substituting these currents into the Hamiltonian Eq. [1.16] we find that the choice 

tanh(277) = (1.18) 

vr 

eliminates cross terms of the form J£(a;) J'j^{x). We now introduce primed Boson fields (j^'Lcix) = 

~7n i'^L] + '^Li) ^^'^ '^'rc associated with the primed charge currents. A factor of \/2 is needed to 



reproduce the correct anomaly: J£(x) = \/2/Tr dx(p'x^^{x) with a similar formula for the right sector. 
The Hamiltonian written in terms of these fields is simply: 

Thus the Bosonic theory remains Gaussian, even for Ac 7^ 0. 

It might be expected that the spin current coupling A^ could be incorporated in a similar fashion. 
However, the spin interaction A^ Ji,g{x) Jro,{x) differs in a fundamental way from the charge 



interaction Ac Jl{x) Jr{x)- Only the JlzJlz part of the interaction has a quadratic representation 
in terms of the Boson field (j)s. The other two components are rather more complicated. Note that 
SU(2) invariance may be employed only when both currents in the bilinear are of the left or right 
type; for example, when we replace Jl0{x)J2^^{x) -^ 3Jl^{x). The non-trivial nature of this term 
is apparent in the Fermion basis: Xg is the only interaction coefficient in Sint that renormalizes. To 
second order in weak-coupling perturbation theory (see below) it flows as: 

Because no fixed points intervene at intermediate coupling in the original lattice Hubbard model 
(which was exactly solved by Lieb and Wu via the Bethe ansatz^) the flow for the continuum 
problem described by Eq. [1.20] is likely to be qualitatively correct at all Ag. The A^ = flxed 
point of Eq. [1.20] (which is stable when approached from the A^ < side) exhibits an enlarged 
SU (2) L <^ SU (2) n symmetry because the left and right spin currents are decoupled. In other words, 
separate SU(2) rotations on the left and right currents leave the fixed point action invariant. The 
strong coupling fixed point, with A^ ^ oo shows no such symmetry; instead a spin gap opens and 
the electrons pair into singlets that require energy to break apart. 

To understand better the role of spin-charge separation in the one-dimensional problem, we 
examine the single-particle Green's function along the fixed line Ag = 0. Following Shankar'*, we 
introduce momentum space regularization by including a convergence factor e~2"lPl dp along with 
the integration measure, and we take the a ^ limit at the end. Now the (equal time) correlation 
functions for the (j)'^^ = -y= [(p'^^ + </'li] ^^Ids are given by: 

GlAx) S< *'io(i)'*'lo(0) - <I>'lM > 
= i-ln " 



47r a + ix 

Grc{x) ^< ^'n,{x)4>'j,M - 4>'l{0) > 

1 a (1-21) 

= -— In — ■ 

47r a — IX 

The spin Bosons (pg exhibit identical correlations. It is easy to restore the time-dependence of these 
correlation functions by using Lorentz invariance (of course with different velocities in the charge 
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and spin sectors). The Bosonization procedure is completed with the observation that Fermion 
operators are equivalent to exponentials of the original unprimed Boson fields: 

^pLa = /- — exp[-j\/47r CpLa] 

tpRa = exp[i\/47r (/)Ra] ■ 

One way to prove relations Eq. [1.22] is by constructing the Fermion currents with the point splitting 
procedure described above. Then spatial derivatives of the Boson fields appear, and the currents 
defined in Eq. [1.10] are obtained. Since (pic = cosh(ry)0'^^ + sinh(ry)0'^^ and (p^ = -K= {4>Lc + 4>Ls) 
etc. we can combine the Bosonization formulas to find the Fermion two-point functions. 

< ^'^lix.t) V'L/3(0,0) > = 2 X — exp{27r[cosh2(r/) G lc[x ,t) ^ ^mk? [t]) Gn^ix.t)]} 

X exp{27rGLs(a;,t)} 

TT 

and 

< Vlf (x, t) V'R;3(0, 0) >= - (x + iv,t)-^l'^{x + ii;,t)-i/2(a;2 + ^Ifiy^ . (1.23) 

Here the anomalous exponent a = sinh (77). The explicit separation of charge and spin refiects 
both the different velocities of the two types of excitations and the remaining interaction Ac in the 
charge sector. 

The path integral picture yields the following free Lagrangian densities^: 

2A 
^ (1.24) 

where the spatial derivatives in these two expressions implicitly include the different velocities 
factors, Vc and Vs- The Bosonization formulas Eq. [1.22] imply that the U{1)l ® U{l)ii symmetry 
operation is effected simply by shifting the left and right Bosons by, in general, different constants. 
Thus, (pc is shifted ( Slc -^ 4>Lc + ^#= and chjic -^ 4>Rc + -iS= ), but chg remains invariant, reflecting 

v27r v27r 



the fact that the symmetry operation acts only on the charge sector. The Lagrangian density Eq. 
[1.24] remains invariant because d^cpc is unaffected by the shift. The Umklapp term, as expected, 
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breaks the symmetry because it is equivalent to adding the term A3 cos(\/8vr <pc) to Lc[4>c] and the 
cosine clearly changes under a shift of (pc by a constant. 

At this point we might question whether the strange form of the propagator somehow eliminates 
the logarithmic divergences that give rise to the RG flow described by Eq. [1.20]. The answer, which 
is no, may be seen easily in the Boson basis where spin-charge decoupling is explicit. There the 
interaction appears as a term proportional to As cos(-v/87r (j)g) added to Ls[(/>s] of Eq. [1.24] which 
drives the RG flow. [Note that the SU{2)L'S>SU{2)[i symmetry at As = now manifests itself in the 
Boson basis as separate shifts by a constant in (pis and (pRs- Actually, full SU(2) symmetry would 
exhibit invariance under three types of rotations corresponding to the three generators of SU(2). 
Abelian Bosonization, however, forces us to choose a spin quantization axis; consequently the 
theory only exhibits explicit symmetry under rotations about that axis. Invariance under rotations 
in the other two directions remains hidden.] In the Fermion basis, the problem is slightly more 
complicated. The logarithmic divergence that drives the flow described by Eq. [1.20] comes from 
a single loop diagram with four external fields that contains two propagators: one for left moving 
fields and one for right movers (see Figure [2]). (Diagrams with two left or two right propagators 
do not yield logarithmic divergences.) The diagram is most easily evaluated in position space. The 
integral to be evaluated is: 

rs poo 1 

I{s) = {\sf / dx / dtix^ + viey^-^'^ix^ + vier^l^ . (1.25) 

Ji J-00 

Since we are performing a perturbative expansion to order (As)^ it is sufficient to set a = (2^)^ + 
O(A^) equal to zero and study whether the logarithmic singularity persists when Vc 7^ Vs- With 
a = the spatial integral yields a hyper geometric function which, when integrated over time, 
indeed produces the desired logarithm. 

Fermi liquid behavior arises only in the special case Vc = Vg and a = as the two-point function 
Eq. [1.23] must have a single simple pole. However, the discontinuity at the Fermi surface remains 
even when Vc 7^ Vg so long as a = 0. To see this, note that the momentum space occupancy is found 
by taking the Fourier transform of the equal time propagator (which we define to be the average 
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of the correlation function evaluated at times t = 0+ and t = 0— in order to specify definite time 
orderings which also cancel out the imaginary component). Use of the Fermion anticommutation 
relations then shows: 

2nL{k) -l = \ j dx e*(^+^/)- {< V'f (x,t = 0+)V'l^(0,0) > + < 4^(x,t = 0-)^i^(0,0) >} . 

(1.26) 

But as t ^ Oib the two velocities disappear from the correlation function, which equals — and 

yields a step function in momentum space only for a = 0. Apparently spin-charge separation and 
the destruction of the Fermi discontinuity are separate issues.^ Both are characteristic properties of 
Luttinger liquids^. In the following we continue to speak of Fermi points and Fermi surfaces. Clearly 
these points or surfaces should now be defined more generally as manifolds of points in momentum 
space at which the zero-temperature occupancy shows non-analytic behavior characterized by an 
exponent a instead of a discontinuity. In particular, near the Fermi momentum the occupancy 
varies as: n{k) ~ n{kf) + C\k — kf\'^°'sgn{k — kf) where C is a non-universal constant that sets the 
momentum scale. It depends on the the momentum-energy cutoff in the interaction Ac. 

The outline for the rest of the paper is as follows. In section (II) we introduce a model with four 
Fermi points in two spatial dimensions that incorporates some key features of the higher dimension 
interacting Fermion problem. A renormalization group solution of the model away from half-filling 
finds a stable fixed line with four-fold U(l) symmetry that naturally generalizes the one-dimensional 
U{1)l "Xi U{l)ji fixed point symmetry. Bosonization of the Fermions at these points suggests a new 
way to look at Fermions in higher dimension and in section (III) we follow this line of thought to 
arrive at full Fermi surface fixed point manifolds with infinite U(l) symmetry. In section (IV) we 
demonstrate that this way of looking at things yields concrete results by calculating non-analytic 
contributions to the specific heat. And in section (V) we rederive the collective mode equation in 
the new framework. Finally, section (VI) contains some discussion and speculations. 
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II. A 2-D FIXED POINT WITH 4-FOLD U(l) SYMMETRY 

We turn now to a simple model for two-dimensional interacting Fermions that illustrates how 
possibly spin-charge separation might occur in spatial dimensions greater than one. The Fermi 
surface of a nearest neighbor tight-binding model on a square lattice serves as inspiration for the 
model. Instead of treating the entire continuum of Fermi points that make up the Fermi line enclos- 
ing the occupied states, we make our problem tractable with the following drastic approximation: 
we treat each of the four sides of the Fermi surface as a single point labeled by ±1 or ±2 (see 
Figure [3]). At half-filling, these points lie respectively at momenta ib(7r/2, vr/2) and ib(7r/2, — '/r/2) 
but away from half-filling the momentum is generally incommensurate with the reciprocal lattice 
vectors. With this simplification, the infinite set of renormalization group equations is reduced to 
a manageable finite set. 

Note that this model differs from models studied earlier by Schulz^ and Dzyaloshinskii^*' that 
focused on the Van Hove singularities at the four corners (0, ibvr) and (ivr, 0) of the half-filled tight 
binding spectrum. Our model also is not equivalent to two coupled parallel Hubbard chains - a 
system studied by Anderson^^, Finkel'stein^^, and others. It differs in that excitations at points 
±2 propagate at right angles in momentum space with respect to excitations at points ±1 whereas 
excitations in the two-chain system always propagate in parallel (or anti-parallel) directions. Con- 
sequently, different marginal interactions and renormalization group equations appear in our model. 
Anderson's analysis^^ of the two-chain problem led him to conclude that weak interchain coupling 
does not change the one-dimensional physics significantly but recent work by Finkel'stein^^ and 
Castellani, Di Castro and Metzner^^ suggests that interchain coupling is relevant and destabilizes 
the Luttinger liquid. 

We choose the model in part because it emphasizes the role that perfect nesting plays in 
destabilizing various fixed points (we elaborate on the nature of these fixed points below). Thus for 
example processes that transfer an electron across the Fermi surface from, say, point 1 to point -1 
receive the same weight as processes that move an electron from point 1 to 2 because the density 
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of states is non-zero only at the four points. Note that the density of states at each point must 
be held constant regardless of the system size. One might be tempted to give each point the same 
weight as the entire line it replaces, but this choice proves uninteresting as quantum fluctuations 
would be suppressed in the thermodynamic limit. Photoemission experiments on the cuprate 
superconductors^^ provide another justification for our model^^. These experiments show that 
hole pockets form around the momentum points ib(7r/2,7r/2) and ib(7r/2, — '7r/2) as the compounds 
are doped away from the insulating limit. Low-energy excitations near these points may play an 
important role in the normal and superconducting phases. 

Obviously van Hove singularities are ignored in our model. They break scale invariance in the 
Gaussian part of the action and therefore cannot be incorporated into the renormalization group 
scheme because the dispersion relation is not linear at those points. While the singularities may 
or may not be important at half-filling, the Umklapp terms drive various instabilities which the 
singularities are unlikely to prevent. In any case, here it is the problem away from half-filling 
that is of most concern to us. The other major limitation of our model - no small angle scattering 
processes - is a severe approximation to physical reality. Small angle scattering is clearly important 
for example in the formation of momentum-space Cooper pairs. We therefore expect, and indeed 
find, unphysical behavior in certain limits (see below). Nevertheless, our calculation suggests that 
stable fixed points describing whole Fermi surfaces, not just points, exist. 

The electron operators Cxa may be written in terms of the continuum fields at the four points. 
We now allow a to take on values 1, ..., n for the SU(n) case. We consider the general problem 
because it enables us to check our calculations more thoroughly for combinatorial errors. It also 
permits us to study the spinless case n = 1. Let u = x + y and v = x — y where x and y are integers 
labeling the coordinates of a point on a lattice with unit lattice constant and place the system in 
periodic box so that 1 < u, v < L. Then at half-filling the lattice electron annihilation operator 
can be rewritten in terms of the continuum fields at the four points as: 

cxa = Me^'^"^' V'ia(t^) + e-^'^"/^ V-ia(«) + e*'^^/^ V2a(^) + e'^^^/^ ^_2^(^)} . (2.1) 
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Away from half-filling we simply replace 7r/2 -^ kf in this formula. With the replacement Eq. [2.1] 
we see that excitations are constrained to move in directions perpendicular to the Fermi "edges". 
In other words, the Fermions cannot move in arbitrary directions, just forward and backwards 
along the lines depicted in Figure [3]. With this definition we can now break up any four Fermi 
interaction into two types of terms: marginal terms that vary smoothly in space and irrelevant 
terms that oscillate rapidly or contain derivatives. Keeping track of just the marginal terms a little 
algebra shows, for example, that at half-filling the Hubbard interaction — (cx Cxa — n/2)^ generates 
the following four Fermi terms: 

(2.2) 

+ (V'l"2V'-la + V'l°lV'-2a + V'l^V'la + V'l"V'2a)^ 
+ (V'i"V'-2a + V'-'lV'2a + V'2°V'-la + 'ip^'ipla)'^ ■ 

More generally, we can write down all possible four Fermi interactions consistent with the symme- 
tries of the SU(n) spin group and the symmetry of the square Fermi surface. Thus the perturbation 



is: 



2a J 



H^nt = ^ j dudv (Ale (Jl J-1 + J2J-2) + Ai, {Jfpf^^ + J!^pf 
+ A2C (Jl + J-l)(J2 + J-2) + A2. (J53 + J'^Ip){JL + J-2a) 

+ A3 Ki^t'i^-la? + {4"i^-2a? + B.C.] 

+ A4 (Vl"^-2aV'f V'-l/? + i^l%i^-la4''i^2(3 + H.C.) (2.3) 

+ A5 (V'l"V'-2aV'l^iV'2/3 + ^!_'^V'-laV'2^V'l/3 + H.C.) 

+ Ae {lpl°'i)-lai^%1p2l3 + i^i°'i^-la1p2^i^-2l3 + H.C.) 

+ A7 [(V'I"^2a)' + ii^l^i^^la)^ + (4>-l«)' + ii^l^i^la)^ + H.C.]}. 

Here we have again introduced the charge and spin currents, now for each of the four points. Note 
that J-ti = J±i{u), J±2 = J±2iv), and likewise for the spin currents. We rescale the interactions 
by a factor of j^ to keep the density of states constant. Not included in the above expression are 
terms like ^ 6vc (Ji)^ and ^ dvg Jia^iB which simply renormalize the charge and spin velocities. 
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For the above on-site Hubbard interaction, the couphng constants take the following values: 

Xic = X2c = {U/n){2-2/n) 

Ais = A2s = -2U/n 

A3 = U/n 

A4 = 2U/n (2.4) 

A5 = 2U/n 

Ae = 2U/n 

Ay = U/n 
but other bare interactions (ie. nearest-neighbor Coulomb or spin exchange) yield other values (see 

below). 

Each of the nine coupling constants in Eq. [2.3] corresponds to a particular process drawn 
in Figure [4]. Unlike the one dimensional case, we see that a number of these marginal processes 
break spin-charge separation. The current-current terms Aic, A2C, Ai^, and A2S respect it, and 
so do the Umklapp terms A3 and A7 (at least for the physical SU(2) problem), but the other 
terms (A5 and two of the Umklapp terms A4 and Ag) are "mixed" processes that scatter both spin 
and charge. Away from half-filling, the term A5 survives and it is this interaction that will draw 
our attention in the following renormalization group analysis. As in one dimension, the model 
possesses global SU(2) spin symmetry. However, it exhibits the four separate U(l) symmetries only 
if A3 = A4 = A5 = Ag = A7 = 0. In other words, only current-current type interactions preserve 
f/(l)i ® U{1)2 f/(l)_i U{1)^2 symmetry. This behavior clearly differs from the one dimensional 
model, which automatically exhibits U{1)l U{1)r symmetry away from half-filling. Again A5 is 
the sole offending term away from half-filling. 

It is a straightforward, though lengthy, exercise to work out the RG flows to second order 
in the coupling strengths by evaluating one loop diagrams with four external Fermi field lines. 
The diagrams are essentially no different from the one we evaluated in one spatial dimension 
(Figure [2]). This is because at the one loop level only diagrams that contain both a 1 propagator 
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(ie. < ^|"(ti,t) ^iq(0, 0) >) and a -1 propagator (or 2 and -2 propagators) yield logarithmically 
divergent contributions. It follows that A2c, ^2s, and A7 do not flow at this order because these 
interactions contain Fermions at points 1 and 2 (or -1 and 2, etc.) so the requisite propagators do 
not appear. Inspection of the diagrams in Figure [4] reveals the physical origin for this decoupling. 
The three interactions A2C, ^2s, and A7 difl"cr from the other terms in that exchange of momentum 
between the two points is forbidden because the two directions are perpendicular. For example, 
consider the momentum-space version of interactions Aic and A2c- Let p denote momentum in the 
u-direction and q be momentum in the v-direction. Then the interactions take the form: 



Ale 
and 



/ ^ Mp) J-i{-p) +1^ J2{q) J^2{-q)} 



\2c L [Ji(p = 0) + J_i(p = 0)][J2(g = 0) + J_2(^ = 0)] . (2.5) 

In contrast to Aic (and Ai^), only the zero- momentum component of the currents couple in the A2C 
(and A2s) terms. 

The remaining six couplings flow as follows [the prime denotes a derivative with respect to 

7rln(s) ]: 

a;, = 8(1 - l/n)(A3)2 + 2(A4)' - (2/n)A4A6 + 2(l/n - 1)(A5)' + 2{\Qf 

A'l, = n{\is? + 4(n - 2)(A3)' - 2A4A6 + 2{\^f + 2n(A6)' 

A'3 = 4AicA3 + 2(n - 1 - 2/n)\is A3 + (A4)^ + 2A4A6 - n{\&f 

(2.6) 
a; = [2Aic - (2/n)Ai, + 4A3]A4 

A'5 = [-2Aie + 2(l + l/n)Ai,]A5 

A'e = [2Aic + 2(n - l/n)Ai, + 4(1 - n)A3]A6 + [4A3 - 2Ai,]A4 . 
We can perform several checks on these equations. First, the equations close: we have not forgotten 

any marginal operators. Second, if we consider only the pair of points 1 and -1, the equations must 

reduce to the known ones in one spatial dimension. By setting A4 = A5 = Ag = one can easily 

check that the remaining equations (A3 is now the Umklapp term mentioned in the previous section) 

do agree with the known results in one dimension. As another test, note that two terms vanish in 
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the physical case n = 2. In particular, A3 no longer couples to Ai^ because the Umklapp term is a 
SU(n) singlet operator only for the special SU(2) case. Finally, what happens when n = 1, the case 
of spinless Fermions? Many terms vanish because they do not exist for spinless Fermions. Thus, 
the spin singlet Umklapp terms A3 = Ay = by the Pauli exclusion principle and of course the spin 
current terms Ai^ = X2S = because there is no spin. Also, As = — A4 because A4 processes can no 
longer be distinguished from Ag ones and A5 = due to internal cancellations present in Eq. [2.3] 
when the spin label is removed. Equations Eq. [2.6] respect this limit. 

The RG flows described by Eq. [2.6] generically flow to large values. The flows are physically 
sensible: at half-filling Umklapp processes generate various instabilities and the system becomes 
gapped in the charge sector when the interactions are repulsive. Attractive interactions, on the 
other hand can lead to superconducting instabilities. The restricted phase space of our model 
obscures the interpretation of these instabilities. For example, the Goldstone mechanism tells us 
that phases of broken SU(n) symmetry exhibit gapless spin excitations. On the other hand, the 
Higgs mechanism suppresses gapless excitations in the charge sector if the U(l) symmetry breaks. 
But our model retains the character of 1+1 dimensional phase space which is not large enough to 
foster broken continuous symmetries. In any case, our failure to treat the van Hove singularities 
and small angle scattering processes is not as important as it might first seem: these processes are 
unlikely to inhibit the formation of instabilities. 

The spinless case n = 1 is an exception. As noted above, we can take Ai^ = A3 = A5 = and 
A4 = —Xq in this case. The flows are described by the Kosterlitz - Thouless equations: 

A'le = 6(A6)2 

(2.7) 
Ag = 2AicA6 . 

Here the fixed line Ae = is stable for Aic < 0. We may interpret the instability at positive Aic as a 
tendency to form a charge-density wave. To see this, note that the next-nearest-neighbor Hubbard 
repulsion ^(rix+x + ^x+y) n^, where rix = CxC^a is the electron occupancy at site x, leads to the 
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following bare continuum couplings at half-filling: 

\ic = \2c = {Ui/n){2 + 2/n) 

Ais = \2s = 2C/i/n 

A3 = -Ui/n 

A4 = (2.8) 

A5 = 

Ae = -2Ui/n 

A7 = . 
(For the case n = 1, A3 should be set equal to zero by the Pauli exclusion principle.) So repulsive 
nearest-neighbor interactions grow, a tendency towards the formation of a charge-density-wave 
sets in and sites on the even sublattice exhibit different charge density than those on the odd 
sublattice. This behavior is consistent with that found by Shankar^ in his functional RG calculation 
for spinless Fermions and therefore lends credibility to our model. On the other hand, Shankar 
finds a superconducting instability for attractive interactions which contrasts with the stability 
shown by our model. We reconcile this difference by noting that the Cooper instability is driven 
by small angle scattering processes that scatter pairs of Fermions of opposite momentum around 
the Fermi surface. Again the phase space for such processes in our model is severely restricted by 
existence of only four Fermi points, and we should therefore not expect momentum-space Cooper 
pairs. Real-space Cooper pairs can arise, however, as we show below. Nevertheless, negative Aic 
corresponds to an attractive interaction and is indicative of a tendency towards the formation of 
superconducting pairs. 

The problem of spinning Fermions away from half-filling is rather more interesting. Setting the 
Umklapp terms A3 = A4 = Ag = Ay = we obtain the reduced set of flow equations: 

A;, = 2(l/n-l)(A5)2 

a;, = n(Ai,)2 + 2(A5)2 (2.9) 

A's = 2[(1 + l/n)Ai, - Ale] A5 . 
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Here we find a stable fixed line defined by Ais = A5 = and Aic > 0. In fact numerical integration 
shows that it attracts flows starting from the repulsive Hubbard couplings given by Eq. [2.4] (see 
Figure [5]). Again, A2c and A2s do not renormalize at second order. Like Aic, these couplings can 
be non-zero along the fixed line. The instability at negative Aic can be interpreted as a tendency to 
form real-space Cooper pairs something like those proposed in the original resonating-valence-bond 
(RVB) theory of Anderson^^ and collaborators. Thus the Fermion spin permits the formation of 
singlet pairs. 

What is the nature of the stable region of the fixed line? First note that it exhibits both 
spin-charge separation and four-fold U(l) symmetry because only current-current type interactions 
remain. The fixed line thus represents a natural generalization of the one-dimensional Luttinger 
liquid and as such motivates the approach to the continuum Fermi surface problem we describe in 
the next section. Let us first look more closely at our solution by transforming to Boson variables. 
At first it seems strange to contemplate Bosonization in spacetime dimensions greater than two. 
But the problem remains essentially one-dimensional; the Fermions at each of the four points are 
restricted to move along lines. In fact it is convenient to introduce complex space-time coordinates 
analogous to those in 1 -|- 1 dimensions. Let u± = uitit and v± = v^it where we remember that the 
velocities, here set equal to one, are in general different in the spin and charge sectors of the theory. 
Obviously the group of conformal transformations in 2 -|- 1 dimensions is finite. Consequently, 
the model does not possess the infinite symmetries of a true 1-1-1 dimensional conformal field 
theory. But it is the current algebra that concerns us most here. It is the essential ingredient that 
permits us to map Fermions onto Bosons and vice versa at each of the Fermi points. Again either 
Abelian or non-Abelian Bosonization works. In this case we choose non-Abelian Bosonization^^ 
for the spin sector by introducing the Wess-Zumino-Witten (WZW) field g^ (the charge sector is 
still described by an Abelian Boson 0). Non-Abelian Bosonization is superior in the sense that 
it explicitly exhibits global SU(n) invariance. With it we can readily classify all SU(n) invariant 
operators. (Abelian Bosonization hides SU(n) invariance because an explicit choice for the spin 
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quantization direction must be made.) 

The Bosonization dictionary again translates currents defined in terms of the Fermi fields into 
Bosonic operators. For example at point 1 we have: 

Jii'^) = \ - (^u 4>i{u) 

V ^ (2.10) 

Here the free fixed point theory with all the A's equal to zero has its charge sector described by a 
free Lagrangian density: 

^c = ^ {{du+cbif + {du_cb.if + {d,^(t>2? + (5._0-2)'} . (2.11a) 

The spin sector consists of a k = 1 WZW action given by Ss[gi,g2,g-i,g~2] = Si[gi +g^i] + S2[g2 + 

g-2] where: 

^ib] = ^ / dudtTr {dug dug^ + dtg dtg^} 

^'! ^7 (2.116) 

+ dudtdz e^^^Tv^g'^d^g g^O.g g^^d^g). 

127r Jy 

Here the second integral, the topological Wess-Zumino term, is defined by extending the domain of 
the g-field from physical two-dimensional (u, t) space-time to a three dimensional volume V with 
space-time coordinates x^ = {u,t,z). The boundary dV of V is taken to be the (u, t) space-time. 
Of course 5*2 [5] is similar in form to Si[g] but the spatial variable v replaces u. The spin sector 
of the free theory displays SU{n)i (gi SU{n)2 ® SU{n)^i (g) SU{n)^2 invariance because the spin 
currents at the four points are decoupled. 

Now the residual fixed line interactions Aic, X2c and A2S can be included by using the Bosoniza- 
tion rules of Eq. [2.10]. In the Boson language, the 4-fold U(l) symmetry operation amounts to a 
shift in each of the charged Boson fields by a constant: (f>i{u,t) -^ (j)i(u,t) + Fi, etc. Since only 
derivatives of the Boson field appear in the action, it continues to manifest four-fold U(l) invariance 
as expected. On the other hand, local SU{n)i (g SU{n)2 (X> SU{n)^i (g) SU{n)^2 invariance is at 
least partly broken by non-zero A2S which couples together the zero-momentum components of the 
spin currents at the four points. 
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Two issues remain to be investigated in our model: First, how are Fermi statistics maintained in 
the Bosonization scheme, now that there are two spatial directions? And second, how do the residual 
fixed point interactions change the character of the Fermi points? We answer these questions by 
constructing a more general framework in the next section. 
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III. LUTHER-HALDANE BOSONIZATION: INFINITE U(l) SYMMETRY 

Encouraged by the renormalization group flows in our model, we now take a leap of faith, 
advocated most recently by Haldane^, and postulate the existence of a similar fixed point, not just 
for the four Fermi points, but rather for a continuum of Fermi points, in other words, a Fermi 
surface. We outline the construction of the currents and the Hamiltonian first and later [in sections 
(IV) and (V)] demonstrate that the framework reproduces well known results. 

To be definite, we study the case of three spatial dimensions; generalizations to other dimen- 
sions are straightforward. We begin with the charge sector and study a smooth Fermi surface 
parameterized by radial vectors S and T that label a fine, locally flat (and rectangular) mesh of 
grid points on the Fermi surface with spacing A « kf between the points. We also place the 
system in a cubic box with sides of length L and use periodic boundary conditions so that the 
momenta are quantized as pm = -^ m where m is a vector with integer components. The most 
general charge Hamiltonian possessing infinite U(l) symmetry may then be written as: 

^- = I E E ^-(^^ T; q) J(S; q) J(T; -q). (3.1) 

S,T q 

The prefactor of 2 compensates for over-counting due to the symmetry of the summand under 
q ^ — q with S ^^ T. The function l^(S,T;q) encapsulates not only the Fermi velocity Vf of 
the non-interacting system but also the residual Fermi-liquid type interactions 1/,' between quasi- 
particles. The current at each point J(S;q) is now defined in momentum space as: 

J(S; q)^Yl ^(S' "^ + q) ^(S; 1^) i^k+q ^"k - 6lo ^k} . (3.2) 

k 

The subtraction of the vacuum charge expectation value rik =< V'k ^/3k > ^^ ^1- [3-2] amounts 
to normal ordering. Our geometric construction of the currents involves tiling the Fermi surface 
with squat rectangular pill boxes at each grid point S. The boxes have dimensions A x A along the 
surface and extend in height it A/2 above and below the Fermi surface (see Figure [6]). The function 
9{S; k) = 1 if k lies inside the box; otherwise it is zero. Thus VfX functions as an ultraviolet energy 
cutoff. 
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The current commutation relations may now be found by direct computation with the use of the 
canonical anticommutation relation { ipq , t/j^p } = S^ <Jq „• First note that currents in different 
patches commute because the Fermion operators that make up the currents are also located in 
different patches. So: 

[J(S; q) , J(T; p)] = <5i,T { ^ ^(^5 ^ + q + p)^(S; k) [0(S; k + q) - ^(S; k + p)](5^+p,onk 

k 

+ J2 ^(S; k + q + p)9{S; k)[^(S; k + q) - 0(S; k + p)] (3.3) 

k 

x(V'£q+pV'ka-<5q+p,0«k)} • 

In one dimension, the index S just labels the left and right Fermi points and the first sum in Eq. 
[3.3] is the usual quantum anomaly. The second sum vanishes in the limit of infinite bandwidth 
(A -^ oo) because in that case 9{q) = 1 except very deep inside or outside the Fermi sea. In this 
limit, matrix elements of the operator {tpj^ , ipka — ^q+pflnk) vanish and we recover the usual 
one-dimensional Kac-Moody algebra. Thus, for the right movers, 

[J{R;q), J{R;p)]=2^Sg+p,o (3.4) 

where the prefactor of 2 comes from the two spins. We recognize this algebra as the momentum 
space version of Eq. [1.7]. 

One might expect that the natural generalization of the currents to two or three spatial di- 
mensions would take the fields to be organized along narrow rays of vanishing thickness radiating 
outward from the center of the Fermi sea. In fact this approach was adopted by Luther in his 
pioneering work on the Bosonization of free Fermions in higher dimension^^ since it reduces the 
higher dimensional problem to a set of simple decoupled 1-1-1 dimensional systems. However it 
is clear that the procedure breaks down when interactions of the Fermi liquid type are included. 
The charge Hamiltonian, Eq. [3.1], couples charge currents in different boxes at positions S and 
T. As the Fermi surface must have non-zero curvature, any wavevector q that lies inside a tube 
at position S, no matter how small, will be accompanied by a wavevector — q that in general does 
not fit inside the tube at a different point T. 
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The problem is avoided with the use of the squat boxes. The price we pay for this new 
geometrical construction is the introduction of several limits which must be carefully taken in order 
to arrive at the correct commutation relations. This delicate series of limits in fact correspond to 
the Fermi liquid theory limits of w ^ and |q| -^ such that -^^ ~^ 0) the so-called w-limit 
which pertains to collective modes rather than quasi-particle scattering^^ . First we require the 
wavevectors q and p in Eq. [3.3] to be small: |q| < jy and |p| < j^ where we take N -^ oo. Thus 
we may think of q as lying within a small sphere inside the squat box - see Figure [6]. The limit 
insures that only the component of q normal to the surface appears in the quantum anomaly. As 
we shall see, it is the normal component that is needed to reproduce the spectrum of low-lying 
excitations in the free Fermion problem. This geometrical result may be obtained by using the fact 
that rik = 2 for momenta k lying deep inside the Fermi sea and zero far outside. With the limits 
|q| < ^ and |p| < -^ we have ^(S; k -|- q -|- p) ~ ^(S; k). The sum in the first term of Eq. [3.3] can 
be done and we have: 

[J(S; q) , J(T; p)] = 2 ,5s,t [ <5^+p,o A^ (^)3 q . ng + error term ]. (3.5a) 

Here hs is the normal vector pointing outward at point S on the Fermi surface and the error term 
is the second sum in Eq. [3.3] which ruins the Kac-Moody algebra because it is not a c-number 
but rather an operator involving the Fermi fields ^p^ and ip. Note that with the above limit on the 
size of q the magnitude of the quantum anomaly is of order -^ (.2^)^- 

Let us estimate now the size of the error term in the commutation relations. It may be estimated 
by replacing the operator (V'k+q+pV'ka — "^q+p o^-k) by 1 — <5q+p o ^^'^ computing the volume of the 
geometrical complement of the intersection of two 6 functions [9{S; k -|- q)—9{S; k + p)] appearing in 
the second sum. Note that the tops and bottoms of the boxes do not contribute because the matrix 
elements are assumed to be zero deep inside or outside the Fermi sea. Only the sides of the pill boxes 
matter. A simple computation then shows that this term is off order Max{|q|, |p|} x A A (^)^ < 
nv^ (^)^. Therefore, the choice of a squat pill box with A = yiVA makes the error term small (of 
order -j^) in comparison to the quantum anomaly. This second limit is equivalent to the "oj-limit" 
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of Fermi liquid theory as ;^ corresponds to A which is now of order v A^ times larger than the 
momentum q. It is satisfying to have this simple geometrical interpretation of the cj-limit of Fermi 
liquid theory. 

The current algebra Eq. [3.5a] can be put into a more familiar form with a Fourier transform 
over the two components of the momentum perpendicular to the Fermi surface normal vector, q^ 
and p_L . Then we obtain: 

[J(S;g|,,xx) , J(T;p|,,yx)] =2 5s,T |^ <5g||+p||,o L'S\^^-y^) . (3.56) 

Here the current algebra is identical to the usual one-dimensional one, Eq. [3.4], except with 
additional labels S, T, xn and yn that "come along for the ride." Thus the well-developed theory 
of one-dimensional Kac-Moody algebra representations^'^ applies equally well to our generalized 
algebra and we can use this machinery to find the spectrum of states. Of course state counting is 
simple in the abelian case but null states appear in representations of the non-abelian Kac-Moody 
algebra. 

What choice of parameters Vc yield the correct spectrum for the charge sector? The non- 
interacting limit is recovered by making the following choice: V^(S,T;q) = ^Vf{S) Q^^ (^It- Here 
the factor of 17 = A^(^)^ cancels the factors of volume appearing in the current algebra Eq. [3.5] 
and the factor of ^ compensates for the 2 due to up and down spins. With the algebra Eq. [3.5] we 
then recover the free dispersion relation uj{S; q) = f/(S) q • ns = vs • q. To scale the interaction 
coefficients properly, we appeal to Fermi liquid theory and note that the current evaluated at zero 
momentum is equal to the occupancy fluctuation operator summed over the interior of the pill box: 

J(S;q = 0)= j; 0(S;k)(5nk (3.6) 

k 

where (5nk = "0^ V'ak— < V^k V'ak >• Therefore, the Fermi liquid interaction is identical to the 
zero-momentum piece of our interaction term: 

^ Y^ /(k, p) 5nk ^«p = ^ E ^^'(^' T; q = 0) J(S; 0) J(T; 0) (3.7) 

k,p S,T 
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if we identify /(ks,PT) = ^'(S,T;0) by assuming that the Fermi hquid interactions depend only 
on the momenta ks and ps at points S and T of the Fermi surface, not the component of the 
momentum perpendicular to the surface. Evidently, a factor of inverse volume (jj) should be 
included in the interaction term, 

ye(S,T;q) = IvfiS) n-' <5i,T + ;^ K'(S,T;q) . (3.8) 

Note that while our theory contains the same Fermi surface interactions as traditional Landau 
Fermi liquids, the form of the interaction is more general than Fermi liquid theory as it depends 
on q, the momentum of the collective excitation. In Fermi liquid theory, the parameters /(k, p) 
do not depend on q and the momenta k and p appearing in /(k, p) are constrained to lie on the 
Fermi surface. A different extension of Fermi liquid theory which relaxes this constraint on k and 
p is described in the next section. Our calculation of non-analytic contributions to the specific heat 
will highlight the difference between these two types of generalizations. 

We see therefore that in higher dimension, as in one dimension, the Bosonization procedure 
puts the free and interacting components of the Hamiltonian on an equal footing, despite the fact 
that the free piece is quadratic in the Fermion operators while the interaction is quartic. This 
simplicity is a result of the current algebra Eq. [3.5] which permits us to express both terms as 
bilinears in the currents. It is however somewhat deceptive because more general quartic terms, for 
instance the A5 interaction in our simplified model of the preceding section, cannot be expressed as 
bilinears in the current operators. Nevertheless, these interactions have a Bosonic representation, 
albeit a more complicated one. We show how to Bosonize general interactions below. 

First we focus on the spin sector. The total Hamiltonian is a given by the sum of the charge 
and spin Hamiltonians. To form the spin Hamiltonian, we define spin currents. In the general 
SU(n) case we have: 

J^-(S;q) ^ J]^(S;k + q) ^(S;k) {^^^^ V^^k - ^ 5^ i^ll^i^,^} ■ (3-9) 



n 
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Like the charge currents, spin currents at different grid points commute, but the non-Abehan 
Kac-Moody algebra governs currents at the same point in the N ^ oo hmit: 



(3.10) 



[J^(S;q) , J]{S;p)] = {62 '^^ " ^ -^^ S]K+p,o ^^^s 

+ 6}jnS;ci + p)-52Jj{S;q + p) 
The physical SU(2) Kac-Moody algebra can be expressed more succinctly as: 

[riS; q) , J\S; p)] = i J^'^l^ q • fig + ie'^'"^ J^(S; q + p) . (3.11) 

The spin Hamiltonian may then be written: 

^^ = ^ j;j;ys(S,T;q) J(S;q). J(T;-q) (3.12) 

S,T q 

where Vg incorporates the Fermi velocity of spin excitations and spin-spin interactions at different 
points on the Fermi surface: 

K(S, T; q) = ^ VfiS) n-^ SIt. + -1 V;'(S, T; q) . (3.13) 

In general it is not possible to exactly diagonalize the Hamiltonian; the non-Abelian nature of the 
algebra precludes this. We encountered this problem in a simpler form in section (II) where Ais, 
the parameter that couples together spin currents on opposing Fermi points, flows by itself (see Eq. 
[2.9]). 

Both the charge and spin currents are invariant under the local U(l) operation which changes 
the phase of all the Fermions inside a given pill box by the same (time-independent) amount F. If 
k lies inside the box centered at grid point S then 

V'a(k) ^ e^^(^) V'a(k) [onlywhen ^(S; k) = 1] (3.14) 

leaves the currents invariant because the ip^ fields, which transform with the opposite phase factor, 
cancel the overall phase change. Thus the Hamiltonian is automatically invariant under the infinite 
U(l) symmetry. The physical meaning of the invariance is clear: the Fermi liquid type interactions 
preserve the Fermion occupancy at each point in momentum space because quasi-particle scattering 
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is suppressed in the N ^ oo cj-limit. The U(l) symmetry just reflects the local conservation of 
particle number. 

Indeed, it is the existence of an infinite number of conservation laws that makes the charge 
sector of the problem solvable. On the other hand, the free Fermion system also exhibits local SU(2) 
symmetry. So it is rather surprising to discover that the spin current interactions in general break 
the infinite local SU(2) invariances down to a single global SU(2) symmetry. The local invariance is 
broken because spin currents at different points on the Fermi surface must rotate together to keep 
the spin Hamiltonian Eq. [3.12] invariant. The special case of purely local current-current coupling, 
14(8, T;q) = 2 Vf{S) Q^^ ^It' ^^ ^^ exception which restores the full local SU(2) invariance. As 
expected, the Hamiltonian is now exactly solvable: the Hamiltonian describes free spin excitations 
propagating at the Fermi velocity. For this special case only the quantum anomaly in Eq. [3.11], 
2(5 f^q • fis, not the ie J'^{S; q + p) term, is important because of the symmetry q -^ — q. The 
factor 2 compensates for the three spin components. The spectrum of states may now be found 
either by simply choosing a spin-quantization axis, or in an SU(2) invariant manner with the use 
of Kac-Moody representation theory.^'' 

Actually, we can find the excitation spectrum when the interactions described by Vg are non- 
singular. In this case, we may treat the spin currents as semi-classical objects: the right-hand side 
of the commutator Eq. [3.11] can be set to zero by rescaling the currents to be of order one. The 
problem resembles the large-spin limit of a quantum magnet since the currents incorporate a sum 
over XQ, » 1 points in momentum space. If we rescale J"(S;q) — > (A$7)^^ J'*(S;q) then the 
rescaled currents obey: 

[J»(S; q) , J\S; p)] = -^ { '^d^'n q • fis + ie^'^ J^(S; q + p)} 

i^^r (3.15) 

as L ^ cx) with A held fixed. The emergence of the classical limit should not be surprising; 
after all. Landau Fermi Liquid Theory is essentially classical in nature. The free dispersion is still 
determined by the quantum anomaly; only the interactions are treated classically by replacing the 
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current operators with their expectation values J"(S; q) -^ {J'^{S; q)) evaluated in the excited state 
of interest. This procedure trivially reproduces the excitation spectrum of Fermi Liquid Theory. 
Note, however, that the classical limit breaks down in the case of singular interactions. For example, 
the first term in Eq. [3.13], 3 ^/(S) ^^^ <^i t' ^^ singular because the factor Q^^ diverges as the 
number of mesh points increases. In this case the quantum anomaly cannot be neglected and in fact 
is needed to reproduce the free dispersion relation. Likewise, any singular spin current interactions 
that couple different patches on the Fermi surface destroy the classical limit: the small anomaly 
cannot be neglected because of the large interaction. In fact this is the generic situation in one 
spatial dimension, where interactions that couple the left and right points are generally of the same 
order as current-current terms that involve only one point. In other words, there is no sense in 
which the Fermi-Liquid type interactions can be smooth when there are just two Fermi points. We 
return to this point in the discussion of section (VI). 

By introducing Boson fields conjugate to the currents, the Fermi fields and interaction terms 
can be Bosonized^^'^. We proceed by analogy to our construction in one-dimension [section (I)] 
and concentrate on spinless Fermions; it is straightforward to include spin via either Abelian or 
non-Abelian Bosonization. We introduce the coarse-grained Boson field (p{S;x.) and the associated 
Boson current in the direction normal to the Fermi surface: 

J(S; x) = Vi^ ns • V(/.(S; x) . (3.16) 

The Boson field is related to the microscopic fields <p{p) by coarse graining over the pill box: 

'^(S; X) = ^ E :rP^ {e'^'-Hp) + e-'^-H-p)} ■ (3.17) 

The microscopic Boson fields satisfy equal-time commutation relations: 

[<A(p) , 0(q)] = 53+q,o • (3.18a) 

Note that the reality of the microscopic fields 0(x) means that i;^(— k) = (/>^(k) and with this in 
mind the commutation relations Eq. [3.18a] take on the more familiar form: 

[cPip] , 0t(q)] = 5l^ . (3.186) 
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Consequently, the coarse-grained fields obey a natural three-dimensional generalization of the one- 
dimensional equal-time commutation relations Eq. [1.8]: 

[0(S; x) , <t>{T; y)] = i^ ^ 4,T ^(^8 • [x - y]) ^^(x^ - y^) (3.19a) 

where again e{x) = 1 for x > 0; otherwise it equals —1. Here x_|_ denotes the two components of x 
that are perpendicular to the surface normal fis. Note that 5^(0) = (2:^)^ which is the area of the 
base of the pill box. Thus, when x_l = y_L we have: 

[0(S; x) , 0(T; y)] = i^ 6^ e(ns • [x - y]) ; (3.1%) 

otherwise the 0(S;x) fields commute. Furthermore, the Boson currents Eq. [3.16] satisfy the same 
U(l) Kac-Moody algebra Eq. [3.5] as the Fermion charge currents Eq. [3.2] (with half the anomaly 
because we have removed the spin index): 

[J(S; q) , J(T; p)] = n 6s,t 5^+p,o Q • ^s (3.20a) 

or in real space, 

[J(S; x) , J(T; y)] = -i n 5s,t L^fts • Vx ^^(x - y) . (3.206) 

Here the Fourier transform of the currents is given by: 

J(S;x)^^ e^'i- J(S;q). (3.21) 

q 

The Hamiltonian Eq. [3.1] then becomes (for spinless Fermions): 

H = 27r I (fx d'y J^ Fe(S, T; x - y) [ng • V(/<(S; x)] [ht • V(/)(T; y)] (3.22) 

■^ S,T 

and Fermi fields are expressed in terms of the Boson fields as: 

^(S; x) = -|L e^i^-- exp{i^ 0(S; x)} (3.23) 

where kg is the Fermi momentum at grid point S. 
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The N-point Fermion correlation functions are reproduced with the use of the Bosonization 
formula Eq. [3.23]. If, for example, we use the operator identity: 

e^ e^ =: e^+^ : exp(ylS + \{A^ + S^)) (3.24) 

then we find that the Fermion two-point function is given by: 

(V^(S;x) V(T;0)) = ^ ^^ e^'^^''^ exp{^ (0(S;x)<A(S; 0) - </>2(S;0))} . (3.25) 

The Boson correlation function can be computed using the relation Eq. [3.17] and the result is: 

Ge(S;z) = (</)(S;x)0(S;O)-02(S;O)) 

il^ , ,ns •x + m, , . , ^ /o o^\ 

= -— ln( : ) ; [x^AI « 1 (3.26) 

47r la 

-^ -oo ; |x_lA| >> 1 . 

Consequently we obtain the correct Fermion correlation function, coarse-grained over the pill box: 

(V't(S;x) V(T;0)) = f- 6l^ T^^^^ i^f 5\^i.) ■ (3-27) 

It should be emphasized that it is the average over the pill box that results in the ^^(xx) term. 

To close the circle (Bosons -^ Fermions -^ Bosons) we form the Fermion current Eq. [3.2]. In 

real space we utilize the point-splitting procedure: 

J(S;x)=:VHS;x)V'(S;x): 

(3.28) 



lime_>o : exp[— i— — — (^(S;x + nse)] exp[i— — — 0(S;x)] : 



27ra n 

then using the operator identity Eq. [3.24] again we obtain: 



17 \/47r Att 

J(S; x) = - — lime^o exp[-z— — {0(S; x + nse) - (^(S; x)}] exp[-2 Gc(S; x)] 

lis • V0(S;x) 

TT 



(3.29) 



which is identical to Eq. [3.16]. A similar calculation shows that the free Fermion Hamiltonian is 
of the same form as the Boson Hamiltonian Eq. [3.22]: 

H'^ = Y. Vf{S) j Sx ^t(S;x)(ns • V)^(S;x) 

^ (3.30) 

= 2^E ^ /"d3x{(ns-V)<^(S;x)}^ 
s •' 
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As it stands, ip fields located in the same patch and at the same perpendicular coordinates anti- 

commute. For example, the Fermion two-point function Eq. [3.27] is odd under the transformation 

X -^ —X followed by complex conjugation which is equivalent to interchanging the creation and 

annihilation operators in a translationally invariant system. However, fields in different patches, 

and fields in the same patch with x_|_ 7^ y_|_, commute: 

{V'(S;x) , ^(S;y)} = ; x_l = y± 

[V^(S;x) , V'(S;y)]=0; x^/y^ (3.31) 

[^(S;x), V^(T;y)]=0; S/T. 
The commutation relations can be transformed into the correct anticommutation relations by in- 
troducing an ordering operator analogous to a Jordan- Wigner transformation^^. Let 0(S) be the 

ordering operator defined by: 

s-i 
0(S) ^ exp{i I Y. "^(T; q = 0)} (3-32) 

T=l 

where the mesh points T have been arranged in consecutive order. To be definite, we could follow 
Luther's prescription and choose the mesh points to begin at some point (the "north pole") on 
the Fermi surface, spiral outwards, and converge at the antipode ("south pole"). It is straightfor- 
ward to check that the combination 'i/^(S;x)0(S) anticommutes with ^(T;y)0(T) when S 7^ T. 
Commuting statistics are still obeyed when the fields are in the same pill box, but this discrepancy 
can be neglected in the continuum limit A ^ 0. Alternatively, a second ordering operator may be 
introduced to implement anticommuting statistics within the pill box. 

Thus we see that charge sector of the semi-classical Landau theory has been replaced by a 
quantum mechanical theory. The Fermi liquid should be thought of as a zero-temperature quantum 
critical Gaussian fixed point with infinite U(l) symmetry and parameters Vc{S, T; q). No longer do 
semi-classical entities like 6ni. appear: these have been replaced by charge current operators that are 
quantized with the Kac-Moody algebras. On the other hand, we have to resort to a semi-classical 
description of the spin sector because the quantum version appears to be intractable. A geometrical 
meaning has been given to the w-limit and a direct connection between the quasiparticle operators 
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and the Boson fields is made via the Bosonization formulas. To exercise the new framework, we 
rederive some well known results in the next two sections. We concentrate on the charge sector to 
illustrate how the quantum theory reproduces these results. 
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IV. THn{T) CONTRIBUTION TO THE SPECIFIC HEAT 

As a concrete application of the proceeding formalism, we calculate the specific heat of an 
interacting Fermi liquid in three spatial dimensions. We obtain a non-analytic T'^ln(r) contribution 
to the specific heat. The existence of such a term is consistent with careful measurements^^ of the 
specific heat in Helium-3. 

We turn off the spin-spin interactions in the following and for simplicity eliminate the spin index. 

As the nonanalytic behavior arises from small momentum processes, it is permissible to treat the 

Fermi surface in a locally flat approximation. Let the surface normal point in the z direction. Then 

the U(l) Kac-Moody algebra can be written as: [J(S;q), J(T;p)] = 2i7 Qz (5| rp (5p_,_ q. These 

commutation relations are equivalent to those obeyed by Bosonic harmonic oscillator creation and 

annihilation operators once we rescale by a factor of the square root of the momentum perpendicular 

to the surface: 

J(S;q) = V-2f^(?. at(S;-q) ; q,<0 

(4.1) 
= ^/2Q q^ a(S; q) ; q^ > 

where [a(S;q), a^(T;p)] = (5| x '^pq- Thus we can find the spectrum by direct diagonalization of 
the Bosonic harmonic oscillator Hamiltonian: 

Hc = ^ Yl V;(S,T;q)(20)g, at(S;q)a(T;q). (4.2) 

S,T q,gz>0 

We again place the system in a box of dimensions L^ and use periodic boundary conditions so 
the momenta are quantized as q^ = -^ m. The Fermi velocity is given in terms of the Fermi 
energy ej by Vf = y^2ef/rn and the number of states at the Fermi surface, A, is given by A = 
Ss Sq ^(S;^) '^nsq,o = ~~;r~ "^^ ~ if ^^ ' Bccausc the pill boxes completely tile the surface we 
have the sum rule: 

EE^(S'q) = ^^- (4-3) 

S q 

The specific heat is computed by using the standard formula for Bosons: 

^- = 4A^E E ^(S;q)-^S)T (4-4) 

^''^^ S q,g.>0 smh (^) 
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where e(S; q) is an eigenvalue of the Hamiltonian Eq. [4.2] which depends on the momentum q as 
well as the index S that labels the vector space of the patches covering the Fermi surface. Let us 
first consider the case of non- interacting Fermions. The eigenenergies of the Bosonized Hamiltonian 
are simply: 

e(S;q)=i'/ q-fis . (4.5) 

The sum over the patch index S and the components of q parallel to the surface just yields the 
number of states at the Fermi surface, A. The sum over the component of q perpendicular to 
the surface can be converted to an integral. Assuming that the temperature is small (so that the 
thermally excited particle-hole pairs lie within the pill box, in other words ksT « VfX) we then 
find: 

= _4_ i2ksTf ^ (4.6) 

AksT^ Vf 27r ^ '^ 

= klT (^) ^ Ll 

^ ^ 27r2 M 

This result is the correct answer for spinless Fermions and of course it should be multiplied by a 
factor of two to account for the spin. It is remarkable that the Boson formula, Eq. [4.4], yields the 
full specific heat. We take it as further evidence that even for spatial dimensions greater than one 
Bosonization reproduces the entire Fermion Hilbert space. 

We now follow Pethick and Carneiro^^ and focus on quasi-particles separated only by a small 
momentum W = ks — kx (ie. |W| << kf) since a consideration of these processes is sufficient 
to demonstrate the existence of non-analytic contributions to the specific heat. Define two vectors 
u = ks + q and u + p = kx — q. The quantity u-p then functions as a small, rotationally- invariant, 
dimensionless expansion parameter. Here the normalized momenta are defined by p = p/p where 
p = |p|. Figure [7] exhibits the geometry of the interaction. The interaction coefficient V^ may be 
expanded in our cylindrical coordinate system. Note that odd powers of u • p do not appear in the 
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expansion because the sum over grid points and momentum eliminates terms odd in p. 

K'(S,T;q)=a + Mu-p)2 + ... 

, 4g2 (4.7) 

= " + ^4^fTw^ + - 

The expansion parameter is controlled in the low-temperature limit which keeps ^2, the particle- 
hole momentum perpendicular to the Fermi surface, small. (Recall that squat pill-boxes force 
|W| >> |q| > \qz\ in the N ^ oo limit.) 

Actually, the interaction differs from the one that Pethick and Carneiro studied: it couples 
particle-hole pairs at points S and T whereas the Pethick-Carneiro interaction couples the occu- 
pancies riu and nu+p. To be precise, the Pethick-Carneiro interaction has the form: 

^ 6 (u • pf 5nu <5nu+p . (4.8) 

u,p 

This interaction cannot be directly expressed in terms of the currents since it involves products of 
distinct occupancies above and below the Fermi surface whereas the current operator evaluated at 
zero momentum, J(S; 0), averages the occupancy operator over the interior of the pill box. There- 
fore a direct connection with the earlier calculation cannot be made. Nevertheless, our purpose 
here is to show how non-analytic contributions to the specific heat arise in the new framework. 
Other terms may make non-analytic contributions; the interaction Eq. [4.7] is the simplest such 
term within our framework. 

To proceed we diagonalize the Hamiltonian with the aid of a Fourier transform from Fermi 
surface patch index S space to X-space. Let 

at(S; q) = -^^ I d^X e^*^'^ at(X; q) (4.9) 



f 



47ri> 

where J d X 1 = —j^ = (number of patches) then 



dX -— —^ -- e(X;q) aT(X;q) a(X;q) . (4.10) 

J J-A/2 ^^ J-A/2 ^^ h ^^ 



Using Eq. [3.8] and Eq. [4.1] we then obtain the eigenenergies: 
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The sum can be converted to a Riemann integral with the substitution A^ X^w ~^ J d^W and 
we find: 

6(X; q) = vjq. " ^ ^ H^^l ^') ■ (4-12) 

In this equation we discard uninteresting terms proportional to b that make additional analytic 
contributions to the specific heat and keep only the logarithmic piece. We treat this term as a 
perturbation and calculate the specific heat to 0(b); then the change in the specific heat 5Cv due 
to the perturbation is: 

A 32vfTT L f°° ^ q'^ \n{q^) 

(4.13) 



^^^ " ^iksT^ (2vr)3 2vr J, ^^^ sinh2(^) 



WALki _,. _. f^ X* 



'\ 4 



^ ;^ TMT) I —TTT^ dx . 
-K-^vj Jo smh (x) 



In the second line we retain only the term containing the T3ln(T) temperature dependence; analytic 
contributions also appear but again these are not interesting. The integral in the second line equals 
7r^/30 so the final result is: -^ « — ^ 6 -^ — THn{T) . Not surprisingly, this result has the same 
form as that found by Pethick and Carneiro^^ as dimensional analysis guarantees this. A direct 
comparison of the coefficient is meaningless however since our interaction is not the same. 



39 



V. COLLECTIVE MODES 

The curvature of the Fermi surface did not play an important role in the calculation of the 
specific heat. In fact we took the Fermi surface to be flat; consequently the Hamiltonian could be 
rewritten as the sum of products of a single creation and a single annihilation operator (see Eq. 
[4.2]). Collective excitations of the Fermi surface, on the other hand, arise from the curvature. It 
is therefore interesting to derive the spectrum of collective modes within the new framework. For a 
curved Fermi surface the Hamiltonian can contain products, for example, of two creation operators, 
so the more general Bogoliubov transformation is required to diagonalize it. 

Again we concentrate on the charged excitations to illustrate the quantum theory. We diago- 
nalize the Hamiltonian Eq. [3.1] by first taking the matrix square root of Vc 

1 1 



K(S,T;q)=^ Fc^(S,U;q) Fc2(U,T;q) , (5.1) 



u 

then we rewrite the Hamiltonian as: 

1 



^ = ^ E E E ^c^(U,S;q) J(S;q)] [J] y.2(u,T;q) J(T; -q)] 

q U S T 



(5.2) 



= ^ E E ^(^' q) '^^(U; -q) . 

U q 

Here we have introduced new charge currents: 

J(U; ci)^Yl ^^ (U, S; q) J(S; q) (5.3) 

s 

and also use the fact the T4 is a real symmetric matrix [ie. T4(S,T;q) = T4(T,S;q)] so therefore 

1 
Vi? is also symmetric. These new currents obey modified Kac-Moody commutation relations: 

[J(S;q) , J(T;p)] = 20 5l^^[V? DV?]{S,T-c^) (5.4) 

where the diagonal matrix 

I)(S,T;q)=5i,Tq-fis (5.5) 

appears naturally in the implicit matrix product on the right hand side of the equation. (The sum 
over the internal indices in Eq. [5.4] has been suppressed for clarity.) We obtain the spectrum by 
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diagonalizing this modified anomaly. Let the eigenvectors u (S;q) and eigenvalues lo (q) of the 

spectrum carry the label A. Suppressing again internal matrix indices and the momentum q we 

have: 

1 1 

2f) Yl [Vc^ DVc^]iS,T) u^(T) =u;^u^(S) . (5.6) 



1 
Upon matrix multiplying both sides of this equation by Vc^ and defining new eigenvectors u (S) = 

1 
^x ^^ (^' '^) ^^(T) we arrive at the collective mode equation (with an implicit sum over repeated 

indices): 

2n Fe(U, S; q)D(S, T; q)ii^(T; q) = lo'^ {cDu^ {V ; q) . (5.7) 

This equation can be rephrased in a more familiar form by writing Vc explicitely as i Vf i7^^ <^i t + 
j-j V^'(S,T;q) and taking the interaction V^ to be independent of q. Dropping the label A and 
the tilde we find the dispersion relation: 

(q • vs - iv) uiS) + ^ E ^ciS, T) (q • vt) u(T) = . (5.8) 

T 

(Recall that vg = w/ng is the Fermi velocity at grid point S.) Now we multiply each term in Eq. 
[5.8] by q- Vs and make another change of variable by redefining u{T) -^ q • vx ""(T) (with no sum 
over T). The result is: 

2Q 

(q • vs - u;) n(S) + q • vs -^ J^ ^/(S, T) n(T) = . (5.9) 

T 

Recognizing that the sum is just a coarse-grained version of the sum over momenta k lying on the 
Fermi surface (FS): 

T keFS 

we see that we have arrived at the collective mode equation. 

Since zero sound excitations involve global distortions of the Fermi surface that slosh Fermions 
back and forth, the curvature of the Fermi surface plays an important role. For example, solving this 
equation for a perfectly spherical Fermi surface with V^{S, T) assumed to be a constant independent 
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of the angle between S and T we find the zero-sound mode: 

cos(9) 
s — cos(p) 

where {9, (/)) are polar coordinates with the polar axis in the q direction and s = ^pr-r- Also implicit 
in Eq. [5.9] is the renormalization of the Fermion mass. Again assuming a spherical Fermi surface, 
we may use Galilean invariance^^ to find the well-known result: 

- = -^ + ^r, (5.12) 

if we identify V^{6) = Yl'iZo fi PiicosO). Finally, collective excitations in the spin sector of our 
theory are given by the corresponding Fermi liquid formula. Apparently our new formulation of 
the Fermion liquid reproduces well known Fermi liquid theory results. 
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VI. DISCUSSION 

In the preceding sections we showed that a simple model of interacting Fermions in two spatial 
dimensions can lead to a fixed point with local U(l) symmetry despite the fact that the bare 
Hamiltonian is only invariant under global U(l) transformations. We also presented a framework 
for the Bosonization of Fermion liquids in higher dimension. We enlarge upon the connection 
between the simplified model of section (II) and the general problem of Bosonization here. First 
it is clear that the model is pathological in the sense that the residual fixed point current-current 
interactions that couple the four Fermi points are singular; for example, Aic (which couples currents 
at opposite points) is typically of the same order as Vc (which couples currents at the same point). 
In the physical case of a continuous Fermi surface in spatial dimensions of two or higher, interactions 
of this type would be equivalent to a current-current coupling of the form: 

VciS, T; q) = ^vcn-^ 61t, + Ai, fl-^<5i,_T (6.1) 

where — T denotes a mesh point directly opposite point T. In the A ^ limit of a finer and finer 
mesh, U~^ oc A~^ -^ co. Thus the second term in Eq. [6.1] amounts to a singular interaction that 
might be expected to destroy Fermi- liquid type behavior^^. 

For singular interactions, however, the connection between the multidimensional Bosonizaton 
and one-dimensional behavior begins to break down for at least two reasons. First, as we noted 
in section (I), Luttinger liquids in one spatial dimension are characterized by the elimination of 
the discontinuity in the Fermion occupancy at the Fermi surface. Consequently, the Fermion 
distribution is smeared out over some energy scale (set by the energy cut-off in the interaction). 
As long as this cut-off is small compared to lattice energy scales, the continuum analysis of section 
(I) holds. In higher dimensions, however, a second energy scale VfX has to be introduced since the 
Kac-Moody algebra is obtained in our construction only in the c<j-limit of A >> A >> |q| where 
A ^ 0. This limit apparently precludes the incorporation of interactions which eliminate the Fermi 
discontinuity. 

43 



We alluded to a second problem with singular spin-spin interactions earlier: the semi-classical 
limit breaks down because the terms on the right-hand side of the spin current commutation 
relations Eq. [3.11] cannot be neglected when interactions diverge in the A ^ limit. A return to 
the original Luther Bosonization prescription using narrow tubes instead of squat pill boxes appears 
to offer a way out of both difficulties. In this case the energy scale VfX need not be introduced; 
we can think of the higher-dimensional problem as a collection of purely one-dimensional theories. 
However, now only singular "tomographic^^" type interactions are permitted: the current in any 
given tube can couple only to itself or to the current in a tube emerging from a point directly 
opposite on the Fermi surface. Actually, the simplified model of section (II) illustrates this problem. 
Interactions A2c and A2s only couple the zero-momentum components of the currents and therefore 
cannot change the excitation spectrum. Nevertheless the special case of tomographic Bosonization 
may exhibit features of interest. 

A separate, but related, problem of interest arises when the velocity of charge excitations differs 
from that of spin excitations. Fermi liquid theory breaks down in this case because the Fermion 
propagator no longer exhibits a simple pole; instead there is a branch cut. Thus the quasi-particle 
weight Z = even though, as mentioned in section (I), a discontinuity in the Fermion occupancy 
remains at the Fermi surface. Since the two velocities are just parameters appearing in Vc{S, T; q) 
and T4(S,T;q) of our theory, we need not restrict ourselves to setting both velocities equal to a 
Fermi velocity Vf as we did in section (III). Anderson has argued that a difference in velocities 
between the two sectors, rather than singular interactions at the antipode, might account for the 
anomolous normal state properties of the copper oxide superconductors^^. It might be interesting 
to explore the consequences of spin-charge separation within the Bosonization framework presented 
here. Other open problems include the incorporation of van Hove singularities and energy gaps 
within the Bosonization framework. It may also be possible to include fluctuations in the Fermi 
surface shape or topology within a renormalization group approach. 

Finally, the Bosonization procedure outlined in this paper may permit the application of semi- 
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classical approximations to the interacting Fermion problem. Semiclassical approximations cannot 
be directly applied to Fermions because the Pauli exclusion principle guarantees that occupation 
numbers are of order one and thus far from the classical limit. Bosonization bypasses this problem 
by replacing the Fermion variables with Bosonic ones. Indeed, semiclassical approximations to 
Bosonized versions of certain one-dimensional problems have been remarkably useful in the past. 
For example, an analysis of quantum spin chains that begins with the weak-coupling Hubbard 
model describing interacting electrons ends up mapping the low-energy theory onto the non-linear 
WZW sigma model. The behavior of this model in the semi-classical limit explains many of the 
known properties of quantum antiferromagnets^^ . 
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FIGURE CAPTIONS 

(1) Typical band structure and the left and right Fermi points of the one-dimensional problem. 
The dashed lines denote filled states. 

(2) A one- loop diagram that gives a second order contribution to the renormalization of A^. All 
diagrams that contribute at this order have a left and a right moving propagator which connect 
the point (x, t) with the point (0, 0). 

(3) The four Fermi points kept in the model with linear dispersion along the lines. The dotted line 
delineates the Fermi surface of the half-filled nearest-neighbor tight binding model. The inset 
shows the (u, v) coordinate system. 

(4) The nine types of marginal interactions: (a) The four current-current interactions that respect 
spin-charge separation, (b) The one non-Umklapp mixed process, (c) The two Umklapp 
processes that transport charge-2 spin singlets and therefore respect spin-charge separation, 
(d) The two other "mixed" Umklapp processes that transport both spin and charge and thus 
break spin-charge separation. (The Umklapp processes only occur at half- filling.) Note that 
these diagrams only depict representative process - the missing diagrams are generated by 
performing the various symmetry operations on the square lattice (reflections that exchange 
points 1 and -1 or 2 and -2 and rotations through 90 degrees). 

(5) Renormalization group flow of the model away from half-filling. The initial coupling is a 
repulsive Hubbard interaction with U /t = 1. The couplings flow toward the fixed line Ais = 
As = and Aic > 0. 

(6) Currents at each grid point S on the Fermi surface are constructed with the use of squat pill 
boxes that tile the surface. The box has dimensions A x A along the surface, height A = -y= 

*= ' *= \/N 

(where N -^ oo) and is bisected through the mid-plane by the Fermi Surface. The function 
0(S;k) = 1 inside the box; otherwise it is zero. The momentum q must be small: |q| < ^. 
Thus, |q| << A << A. 
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(7) Geometry of the Fermi liquid interactions that lead to non-analytic contributions to the specific 
heat. Two squat boxes lie on the locally flat Fermi surface at grid points S and T (see text). 
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